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^ ' Abstract. Deeply inside dense molecular clouds and protostellar disks, the interstellar ices are protected from stellar energetic 

C/3 , UV photons. However, X-rays and energetic cosmic rays can penetrate inside these regions triggering chemical reactions, 
molecular dissociation and evaporation processes. We present experimental studies on the interaction of heavy, highly charged 
and energetic ions (46 MeV ^**Ni'^+) with ammonia-containing ices H20:NH3 (1:0.5) and HtOiNHsiCO (1:0.6:0.4) in an 

, attempt to simulate the physical chemistry induced by heavy ion cosmic rays inside dense astrophysical environments. The 

^ ' measurements were performed inside a high vacuum chamber coupled to the IRRSUD (IR radiation SUD) beamline at the heavy 

l/^ I ion accelerator GANIL (Grand Accelerateur National d'lons Lourds) in Caen, France. The gas samples were deposited onto 

0^ ■ a polished Csl substrate previously cooled to 13 K. In-situ analysis is performed by a Fourier transform infrared spectrometer 

' (FTIR) at different fluences. The averaged values for the dissociation cross section of water, ammonia and carbon monoxide due 

, to heavy cosmic ray ion analogs are ~ 2x10"'^, 1.4x10"'^ and 1.9x10"'^ cm-, respectively. In the presence of a typical heavy 

■ cosmic ray field, the estimated half life for the studied species is 2-3 xlO* years. The ice compaction (micropore collapse) due 
T-H ' to heavy cosmic rays seems to be at least 3 orders of magnitude higher than the one promoted by (0.8 MeV) protons . The 
0^ . infrared spectra of the irradiated ice samples present lines of several new species including HNCO, N2O, OCN" and NH^ . In 

' the case of the irradiated H20:NH3:CO ice, the infrared spectrum at room temperature reveals five bands that were tentatively 

^ I assigned to vibration modes of the zwitterionic glycine (NH3CH2COO"). 
•i-H . 

' Key words, astrochemistry - methods: laboratory - ISM: molecules - Cosmic-rays: ISM - molecular data - molecular pro- 

^ ' cesses 
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1. Introduction sis of such ices can create complex organic compounds, and 

even pre-biotic molecules like amino acids and nucleobases 

The birthplace of stars is the densest and coldest place of the in- ^ Bernstein et al. 2002; Munoz Caro et al. 2002; Kobayashi 

terstellar medium (ISM) called dense molecular clouds. These 20O8) 
regions have typical gas densities of 10^-10^ atoms cm"^ and 

temperatures of the order of 10-50 K. Due to the low temper- The observation of molecules in the gas phase deeply inside 
ature, dust particles (mainly silicates and carbonaceous com- these cold and dense regions, where the gas sticking efficiency 
pounds like amorphous C and SiC) can accrete molecules from on grains is close to unity, suggests that they are indeed energet- 
the gas phase and get coated with an ice mantle. The interstellar ically active regions (e.g. Ehrenfreund & Charnley 2000). Only 
ice mantles (or simply interstellar ices) are composed primarily a negligible amount of stellar UV photons reaches the inner 
of amorphous H2O but usually also contain a variety of other parts of dense regions (Tielens & Hagen 1982; d'Hendecourt 
simple molecules such as CO2, CO, CH3OH and NH3 (e.g. et al. 1985). Other mechanisms have therefore been proposed 
Ehrenfreund & Charnley 2000; Ehrenfreund & Shuttle 2000; to explain the presence of free molecules such as cosmic-ray- 
Gibb et al. 2001; Boogert et al. 2004). Laboratory studies and induced UV photons (Prasad & Tarafdar 1983) and direct cos- 
astronomical observations indicate that photolysis and radioly- mic ray interaction with ice mantles (e.g Shen et al. 2004). Both 

processes lead to molecular desorption from the surface. As 

Send of/print requests to: S. Pilling, pointed out by Shen et al. (2004), in the case of cosmic ray par- 
e-mail: sergiopilling@yahoo.com.br tides, the sputtering promoted by direct impact and the whole 
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grain heating (classical sublimation) due to energy deposition 
in the bulk are also two efficient processes to release molecules 
from frozen surfaces to the gas phase. 

The first detection of ammonia ices has occurred around the 
young stellar source (YSO) NGC 7538 IRS9 (Lacy et al. 1998). 
It has also been observed around the massive protostar GCS3 
(Chiar et al. 2000) and several other massive YSO like W33A 
(Gibb et al. 2000; 2001). Ammonia has also been detected ex- 
tensively inside solar system in cometary coma and on several 
moons (e.g. Kawakita et al. 2006; Bird et al. 1999; Moore et al. 
2007 and references therein). 

One of the most reactive and important molecules observed 
from photolysis and radiolysis experiments of ammonia- 
containing ices, in which there is a carbon source, is the cyanate 
ion, OCN" (e.g. Demyk et al. 1998; Hudson & Moore 2000; 
van Broekhuizen et al. 2005 and references therein). Its in- 
frared band around 2165 cm"' was observed in several proto- 
stellar sources (e.g. Lacy et al. 1984; Gibb et al. 2000; Whittet 
et al. 2001) indicating regions were UV photochemistry and 
ion/electron bombardments have played an important role in 
the chemical alteration of grain mantles. 

Although the flux of heavy ions (e.g. Fe, Ni, Si, Mg, ...) 
is about 3-4 orders of magnitude lower than that of protons 
(Roberts et al. 2007; Mennella et al. 2003), their effects play 
an important role on interstellar grains since they can deposit 
about 100 times more energy than the light ions (He^ and pro- 
tons) inside the grains. Brown et al. (1984) and Fama et al. 
(2008) using high energy He^ and protons have found that 
the sputtering yields scales with the square of the electronic 
stopping power. We have recently proved that this square law 
extends its validity up to values of electronic stopping, such 
as those induced by swift heavy ions (Seperuelo Duarte et al. 
2009a,b). Consequently, amount of species released per impact 
to gas phase due to heavy ions could be 4-5 orders of magnitude 
higher than those for protons. 

In order to characterize the effects of heavy ions on dif- 
ferent interstellar ice analogs, it is necessary to measure their 
sputtering yields and chemical reaction cross sections and to 
compare them to those promoted by He^ and protons. Studies 
on the effect of heavy projectiles in astrophysical ice analogs 
are scarce. Most of experiments have been performed in the 
keV or hundred of keV range, mainly with protons. He and Ar 
ions (e.g. Loeffier et al. 2006; Gomis et al. 2004a,b). Due to 
the low kinetic energy of these ions, the structural and chemi- 
cal changes in the frozen ices are mainly governed by nuclear 
energy loss in elastic ion-target atom collisions (e.g. Brown et 
al. 1982; 1984). 

In this work, we present infrared measurements of 
ammonia-containing astrophysical ice analogs (H20:NH3 
(1:0.5) and HjOiNHjiCO (1:0.6:0.4)) iiTadiated by 46 MeV 
58]yji3+ simulate the modification induced by heavy cosmic 
rays inside dense astrophysical environments. At such high ion 
velocity, the energy deposition is mainly due to inelastic elec- 
tronic interactions with target electrons (electronic stopping 
power regime). 



to pump 




Sample on I ■ I Vacuum 

Csl window jon beam chamber 



Fig. 1. Schematic diagram of the experimental set-up. The ion 
beam impinges perpendicularly on the thin ice film deposited 
on a Csl crystal. 

2. Experimental 

We have used the facilities at the heavy ion accelerator GANIL 
(Grand Accelerateur National dTons Lourds) in Caen, France. 
A 46 MeV ^**Ni'^^ ion projectiles impinge perpendicularly on 
the ice target. In this paper, results for two ammonia-containing 
ices H20:NH3:CO (1:0.5:0.4) and H20:NH3 (1:0.5) are pre- 
sented. The incoming charge state 13-1- corresponds approxi- 
mately to the equilibrium charge state after several collisions 
of a 46 MeV Ni atoms (independent of the initial charge state) 
with matter (e.g. Nastasi et al. 1996). 

In-situ Fourier transformed infrared (FTIR) spectra were 
recorded for ices irradiated at different fluences, up to 2 x lO'^ 
ions cm"- using a Nicolet FTIR spectrometer (Magna 550) 
from 4000 to 650 cm"' with 1 cm"' resolution. The ion flux 
was 2 X 10^ cm"^ s"'. A background allowing the absorbance 
measurements is collected before gas deposition. The total ir- 
radiation time was about 3 hours for each sample. For further 
experimental details see Seperuelo Duarte et al. (2009a). 

The sample-cryostat system can be turned over 180° and 
fixed in three different positions to allow: i) gas deposition, ii) 
FTIR measurement and iii) perpendicular irradiation as shown 
in Fig [T] The thin ice film was prepared by condensation of 
gases (purity better than 99%) onto a Csl substrate attached 
to a closed-cycle helium cryostat, cooled to 12-13 K. A very 
small fraction of the mixture CO2, CO'^O and C'**02 existed 
as contaminants in the pre-chamber. During the experiment the 
chamber pressure was around 2 x 10"** mbar 

Following d'Hendecourt & Allamandola (1986) the molec- 
ular column density of a sample was determined from the rela- 
tion between optical depth Ty - ln(Io/I) and the band strength, 
A (cm molec"'), of the sample respective vibrational mode. In 
this expression / and Iq is the intensity of light at a specific 
frequency before and after passing through a sample, respec- 
tively. Since the absorbance measured by FTIR spectrometer 
was Absy = log(/o//) the molecular column density of ice sam- 
ples was given by: 

1 r 23 r 

N - — \ Tydv - — Absydv [molec cm"^] (1) 
A J A J 

where Absy = Inik/I)/ In(lO) = Tv/2.3. 
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Table 1. Infrared absorption coefficients (band strengths) used 
in the column density calculations for the observed molecules. 



Frequency 


Wavelength 


Assignment 


Band strength (A) 


Ref. 


(cm-i) 


(jim) 




(cm molec"') 




2342 


4.27 


CO" (V3) 


7.6x10-" 


[1] 


~ 2234 


4.48 


NjO (vj) 


5.2x10-" 


[2] 


-2165 


4.62 


OCN- (yj) 


4x10-" 


[3] 


2139 


4.67 


CO (vi) 


1.1x10-" 


[4] 


~ 1100 


9.09 


NH3 (yj) 


1.2x10-" 


[5] 


~ 800 


12.5 


H2O (yz.) 


2.8x10-" 


[4] 


" C0"*0 and C"*02 at 2325 and 2309 cm"', respectively. 



[1] Gerakines et al. 1995; [2] Wang et al. 2001; [3] d'Hendecourt & 
AUamandola 1986; [4] Gibb et al. 2004; [5] Kerkhof et al. 1999. 

From these measurements and assuming an average density 
for the ice samples of about 1 g/cnv' (water ice), the thickness 
and the deposition rate are determined. For the H20:NH3:CO 
ice (1:0.6:0.4) mixture the thickness was of about 1.7 jum and 
the deposition rate of ~ 13 jum/h. For the H20:NH3 ice (1:0.5) 
the thickness was ~ 1 .4 /im and the deposition rate was about 
10.5 /im/h. 

The analyzed ice layers were thin enough: i) to avoid sat- 
uration of the FTIR signal in transmission mode and ii) to be 
fully crossed by ion beam with approximately the same veloc- 
ity. This latter point is important since a relatively small total 
kinetic energy loss of the projectile in the film guarantees that 
the studied cross sections remain constant. 

The vibrational band positions and their infrared absorption 
coefficients (band strengths) for used in this work are given in 
Tabled 

3. Results 

Figure |2^-b presents the infrared spectra of H20:NH3:CO ice 
(1:0.6:0.4) before (highest curve) and after different irradiation 
fluences. Each spectrum has an offset of 0.05 for better visu- 
alization. The narrow peak at 2100 cm-' is the CO stretch- 
ing mode (vi). The broad structure from 3100 to 3500 cm-' 
presents a combination of between vibration modes of water 
(vi) and ammonia (vi ). The band at 1600 cm-' is composed by 
two lines corresponding to the water V2 vibration mode (1650 
cm'') and ammonia V4 vibration mode (1630 cm-'). The fea- 
ture around 1100 cm-' is the umbrella vibration (V2) mode of 
ammonia and the one at 800 cm-' is the libration mode (v/.) of 
water molecules. The OH dangling bond (OH db) line at about 
3650 cm-' is also observed. This band indicates a high degree 
of porosity (Palumbo 2006) in this ice mixture and will be dis- 
cussed later The region between 2400 to 1200 cm-' is shown 
in detail in Figure|2j3. We can see the (V3) band of newly formed 
cyanate ion OCN- at 2160 cm-' and NH^ around 1500 cm-'. 
Increasing of CO2 abundance as a function of fluence is mainly 
due to following processes (CR symbolized cosmic rays): 

CO C + O (2) 

CO + O CO2 (3) 



H^OiNH^iCO ice (1:0.6:0.4) 




4000 3500 3000 2500 2000 1500 1000 

Wavenumber {cm"^) 




2 4 6 8 10 12 14 16 

Fluence (10'^ ionscm"^) 



Fig. 2. (a) Infrared spectra of H20:NH3:CO ice (1:0.6:0.4) be- 
fore (top dark line) and after different irradiation fluences. (b) 
Expanded view from 2400 to 1200 cm-'. Each spectrum has 
an offset of 0.05 for better visualization, (c) Molecular col- 
umn density derived from the infrared spectra during the ex- 
periment. 

The variation of the column density of the most abundant 
molecules observed during the iiTadiation of H20:NH3:CO ice 
(1:0.6:0.4) by 46 MeV Ni ions is shown in Fig.|2}; as a function 
of fluence. The column density of water was determined taking 
into an account the integrated value of IR absorbance between 
1025 cm-' and 800 cm-' (approximately the half left part the 
IR libration band) and multiplying it by a factor 2. Since this 
procedure was done for all IR spectra, the relative error is ex- 
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pected to be small. In the case of CO2, the column density is 
the sum of all its isotopic contributions. For both CO'^O and 
molecules, the same value of band strength determined 
for the V3 vibration mode CO2 at 2342 cm"' (Gerakines et al. 
1995) was adopted. 

The ammonia and CO abundances seem to present a simi- 
lar dissociation rate reaching the half of initial values at a flu- 
ence of about 4 xlO'^ ions cm"^. The column density of wa- 
ter is constant, about 2 xlO'^ molecules cm"^, decreasing very 
slowly as the fluence increases. This is attributed to a persistent 
deposition of water from the residual gas. This effect was not 
observed in our previous experiments on pure CO and CO2 ices 
(Seperuelo Duarte et al. 2009a,b) and is then apparently related 
to the deposition of water containing ices.) 

The OCN" abundance increases very fast, reaching a max- 
imum at a fluence between 1.2 and 1.3 xlO'-' ions cm"^. The 
N2O production reaches a maximum at around 5 xlO'^ ions 
cm"^ and remains constant until the end of irradiation (up to a 
fluence of 2 xlO''' ions cm"^). As pointed out by Jamieson et 
al. (2005), N2O was observed in the radiolysis of N2-CO2 ices. 
They have proposed a pathway reaction between N2 and the 
oxygen atom coming from the dissociation of CO2 by 5 keV 
electrons. N2O was also reported in both N2-HCO and N24-CO2 
ice irradiations by 0.8 MeV protons from a Van de Graff ac- 
celerator performed by Hudson & Moore (2002). Recently, we 
performed additional experiments on N2-CO (1:1) and NH3- 
CO (1:1) ices irradiated by 537 MeV Ni projectiles. N2O 
species have been observed for the N2-CO system, confirming 
the results reported by Jamieson et al. (2005) and by Hudson 
& Moore (2002), but not for NH3-CO. The conclusion is that 
N2 molecules are not easily formed from NH3 and that they are 
required as an intermediate step to synthesize N2O. 

4. Discussion 

4.1. Dissociation cross section 

In the present case of fast ions (projectile velocity > 0.2 
cm ns"'), most of the deposited energy leads to excita- 
tion/ionization of target electrons. In turn, the electrons liber- 
ated directly from the inner part of the ion track, about 0.3 nm 
of diameter for 46 MeV Ni ions (Iza et al. 2006), transfer most 
of their energy to the surrounding condensed molecules (~ 3 
nm). The re-neutralization of the track proceeds concomitantly 
with the local temperature rise, leading to an eventual sublima- 
tion. At the ice surface, the energy delivered by these heavy 
ions allows molecules to be removed by sputtering, causing the 

Table 2. Stopping power and penetration depth values for 
equivelocity 46 MeV ^^Ni and 0.8 MeV protons in water ice 
calculated by SRIM code. 



ion 


Stopping Power 


Penetration depth 




(keV/yum) 






electronic nuclear 




0.8 MeV p 


29 0.02 


111 


46 MeV Ni 


4900 15 


19.5 



formation of craters of a few nm deep (Schmidt et al. 1991). 
The sputtering yield, Y, scales with the square of the electronic 
stopping power inside the target, S e -dEjAx (e.g. Brown et al. 
1984; Fama et al. 2008, Seperuelo Duarte et al. 2009a,b). 

The electronic and nuclear stopping power values for 46 
MeV ^^Ni and its equivelocity protons (0.8 MeV protons) in 
water ice (p = 1 g cm"-'), calculated by SRIIvfl- Stopping and 
Ranges of Ions in Matter code (Ziegler & Biersack, 2006), are 
given in Table |2] Their penetration depths in the water ice are 
also shown. The electronic stopping power is at least 2 orders 
of magnitude higher than the nuclear one; the total stopping 
power for 46 MeV Ni atoms is about 2 orders of magnitude 
higher than that of equivelocity protons. 

The variation in the molecular abundance due to the in- 
coming radiation can be attributed to several processes: i) the 
molecular dissociation quantified by the dissociation cross sec- 
tion, (Td (cm^); ii) the sputtering yield, Y (molecules desorbed 
per ion impact); iii) the formation of molecules, crj (cm^), from 
species already present as well as from radicals and ionic frag- 
ments; iv) the molecular layering, L (molecules deposited per 
ion impact), due to deposition of the residual gas or the even- 
tual recondensation of sputtered molecules. The layering yield, 
L, is expressed in terms of molecules deposited by ion impact; 
it is obtained by dividing the layering rate (molec s"') by the 
ion projectile flux (2 x 10'' ions s"'). 

Taking into account the two "formation" processes (chem- 
ical reaction and layering), as well as the two "disappearance" 
processes (molecular dissociation and sputtering, the column 
density rate for each molecular species is given by: 

^ = ^ cTfjjNj + Li - cTdjNi - YiQiiF) (4) 

where Yjj '^fjjNj represents the total molecular production rate 
of the / species directly from the j species, F, is the sputtering 
yield of a pure ice formed by / species, and Q,(F) is the relative 
area occupied by the / species on the ice surface. 

Assuming: i) that only / - 1 species (e.g. water) of the resid- 
ual gas are able to condense on the sample (i.e., L,- = Lidn, 
where d/j is the Dirac function), and ii) that the analyzed molec- 
ular species cannot react in a one-step process to form another 
species originally present in the ice (i.e., o-fjj x 0), this equa- 
tion system becomes: 

dN- 

= -(rd,iNi - (F,Q,(F) - Li5,i) (5) 

dF 

As shown in Fig. [3] and Fig. |4^ (hachured region), there is 
a non-exponential decreasing behavior at the beginning of the 
irradiation, which boundary is defined hy F < F,„i„ =s 5 x lO" 
ions cm"^. This suggests that ice compaction is occurring, as it 
will discussed further. 

Another phenomenon that can also occur is layering of 
species 1 over the ice surface (Q,(F) — > dn), preventing pro- 
gressively sputtering of / 1 species. At the end of these pro- 
cesses, the equation system (|5]l is decoupled and can be solved 
analytically: 

A^,- = A^o; e\p(-o-djF) , for / 9^ 1 species (6) 
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Table 3. Dissociation cross sections of the studied molecular species for radiolysis of ammonia-containing ices by 46 MeV Ni. 



Species Mixture o\i T/Z Y U A'o Model 
(HzOiNHjiCO) (10-" cm") (10" molec cm-^) (10'' molec ion"') (10-* molec ion-') (10" molec cm-^) 

hIo (1:0.5:0) ^ 23 P 38 29 1 

(1:0.6:0.4) ~2 19 1" 35 24 2 

NH3 (1:0.5:0) 1.3 NA' 0* 2.0 3 

(1:0.6:0.4) 1.4 NA 0* 1.7 4 

CO (1:0.6:0.4) 1.9 NA 0* 1.0 5a 

(1:0.6:0.4) 1.9 NA 0* I'' 1.0 5b 

(1:0.6:0.4) 1.9 NA V 1.0 5c 



" Taken from Brown et al. (1984). * No sputtering. Assuming that the water layering is thick enough to fully cover the NH3 or CO molecules 
on the surface. Assuming that there is no water layering and considering a sputtering of 1x10* molec ion"' (the same adopted for water 
molecules). Assuming an extra source of CO for Layering. ' NA=Not applied. ^ For H2O: derived from L = N^crj + Y. For NHj and CO: 
assuming no layering (except model 5b). 



o H^OiNHjiCO ice (1:0.6:0.4) * H^OiNHj ice (1 :0.5) 




Fiuence (10^^ ions cm ^) 



Fig. 3. Variation of the column density of water, ammonia and 
CO as a function of fluence. Solid lines are the fittings using 
Eq. |6] and Eq. |7] Triangles and circles correspond to different 
ice samples. Model parameters are listed in Table[3] 

and 

A^i = {Noi - A^ci) expi-a^jF) + N^^i (7) 

where A^oo,i = (ii - Yi)/(Td,i is the asymptotic value of column 
density at higher fluences due to the presence of layering. A^, 
and A'o, are the column density of species / at a given fluence 
and at the beginning of the experiment, respectively. Eq. [Tide- 
scribes the condensation of / = 1 species during irradiation. 
This equation is particulary useful for species which adsorbs 
on the walls of vacuum chamber as the case of water, generat- 
ing a continuous flow towards the cold target, creating on it a 
superficial layer during the experiment. 

Fig. [3] presents the best fittings of the infrared data of ices 
H20:NH3:CO and H20:NH3, employing equations |6] and [T] for 
describing the column density evolution with the fluence. The 
condition cry iy ^ applied for water means that the pro- 
duction rate of water molecules from ammonia and carbon 
monoxide is negligible compared to the layering rate. Since 



F„,i„ ^ 5x10" ions cm"^ for the current data, the fittings cover 
only the F = 1 - 12 x 10'^ ions cm"^ range to avoid both the 
mentioned transient effects and eventual large inhomogeneities 
in the ice induced at too large fluences. The observed decreas- 
ing of A', for CO and NH3 ices is indeed exponential, leading 
to the CO dissociation cross section the value of 1.9 x 10"'^ 
cm^, in agreement with the value obtained for pure CO ice irra- 
diated by the same projectiles (Seperuelo Duarte et al. 2009b). 
The values for NH3 are 1.3 and 1.4 xlO"'^ cm^, obtained with 
the H20:NH3 and H20:NH3:CO ice, respectively. The fitting 
parameters (dissociation cross section, sputtering, layering and 
the initial relative molecular abundance of each species in the 
ices) are listed in Table [3] The individual effects of layering or 
sputtering on the column density (e.g. for CO species) coiTe- 
spond to the models 5b and 5c, respectively. 

The H2O column density data, presented in Fig. |2j;, level 
off around 2 x 10'^ molec cm"^ for the two samples. The aver- 
age value for the dissociation cross section of water, employing 
Eq|6] is cTj 1 ~ 2 X 10"'^^ cm^. The sputtering yield for water 
measured by Brown et al. (1984) was extrapolated for the 46 
MeV Ni ions impact by using the stopping power displayed 
on Table |2] The obtained value is Fi = 1x10"* molecules 
per impact; the estimated average water layering of both ex- 
periments is 4 X 10^ molec ion ', obtained from the relation 
L = NooCTd + Y. Due to the large value of water layering, the 
NH3 and CO species were recovered by a H2O film during irra- 
diation and their sputtering yields were considered negligible; 
therefore, the data were adjusted directly by Eq|6l 

4.2. Ice compaction 

A high density amorphous form of water ice is found when 
water vapor is deposited on very cold (T < 30 K) surfaces, at 
a rate < 100 yum/h (Jenniskens, Blake & Kouchi 1998). Bofli 
of our experiments are within this regime, nevertheless the ice 
mixture containing CO was found to be less compact than the 
H20:NH3 ice. This conclusion is consistent with the observa- 
tion of large changes in the column density of the adsorbed 
species (Fig. |4^) and is supported by the presence of the OH 
dangling bond (OH db) line around 3650 cm ' (Fig.|4}5), which 
is attributed to water molecules in the surfaces of micropores 



6 



Pilling et al.: Radiolysis of ammonia-containing ices by heavy and energetic cosmic rays 



o HjOiNHjiCO ice (1:0.6:0.4) » HjO:NH3 ice (1:0.5) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Fluence (10'^ ions cm ^) 
HjO:NH3:CO ice (1:0.6:0.4) 
1.0 Fp — • — ' — ■ ' ' ' ' ' ' ' ' ' ' q 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Fluence (10'^ ions cm"^) 



Fig. 4. (a) Expanded view of Fig. [3] at the beginning of irradi- 
ation. Note the drastic changes in column density of water and 
ammonia up to fluence of 5x10" ions cm"^. Hachured strip 
indicates the regions where the compaction is occurring, (b) 
OH dangling bond feature area (cm ') in the H20:NH3:CO ice 
(1:0.6:0.4) as a function of fluence. Figure inset is the OH db 
feature from non irradiated ice and at fluences 0.1, 0.5, 1, 2.5 
and 5x10" ions cm"^. 



(Rowland & Devhn 1991; Rowland et al. 1991; Palumbo 2006 
and references therein). Since both ices were deposited roughly 
at the same rate (~ 12 fim/h) and at the same temperature 
(~ 13 K), the porosity enhancement must be attributed to the 
CO molecules, which diffuse very efficiently into the ice (e.g. 
Collings et al. 2004). Some cavities are eventually created dur- 
ing the ice formation reducing its compaction. 

Figure |4^ presents an expanded view of Figure [3] corre- 
sponding to fluences up to 3.5x10'^ ions cm"^, i.e. at the begin- 
ning of the irradiation. The hatched region (F < F„,i„) indicates 
the fluences for which the column densities present a complex 
behavior Such rapid variations, observed for both water and 
ammonia species, are enhanced in the case of H20:NH3:CO 
mixture but are not seen in the CO column density. This fact 
may be associated with changes in the band strength of vibra- 
tion modes of hydrogen-bound molecules produced by com- 
paction efi^ects at the beginning of the irradiation. 



A direct evidence of the compaction (micropore collapse) 
due to the impinging heavy ions is shown in Figure |4j3, where 
the intensity of the OH db band as function of fluence is exam- 
ined. The experimental data (band area, y) have been fitted by 
the exponential curve y oc e where F is the ion fluence (ions 
cm"^) and cr = 1 x 10"" cm^ is the disappearance cross section. 
The estimated error is within 30%. The figure inset presents the 
OH db feature at 3650 cm ' during the beginning of the irra- 
diation. The traces correspond to fluences (non-irradiated), 
0.1, 0.5, 1, 2.5 and 5x10" ions cm"^. The peculiar changes in 
the column density, at the beginning of the irradiation, and the 
fact that the OH db feature seems to disappear approximately 
at the same fluence, around 5x10" ions cm"^, indicate that at 
this point the ice is fully compacted by the Ni ion beam. 

A detailed investigation on the compaction elTects due to 
ion bombardment was performed by Palumbo (2006) using 200 
keV protons in water-rich ices at 15 K. The intensity of the OH 
db line decreases after ion irradiation and the amount of ab- 
sorbed carbon monoxide (from an upper layer of CO ice) also 
decreases as the fluence of impinging ions increases. These 
results indicate that the porosity of amorphous water ice de- 
creases after ion irradiation. Palumbo (2006) has quantified the 
compaction degree of water ice irradiated by 200 keV protons 
by fitting an exponential function to the integrated area of the 
OH db peak at different fluences obtaining a "disappearing" 
cross section of 4.13 xlO cm^. Comparing the cross section 
of both experiments, it is observed that compaction effects due 
to heavy and energetic ions (46 MeV Ni ions) are 3 orders of 
magnitude higher than the one promoted by (0.8 MeV) protons. 

4.3. Synthesis of molecules 
4.3.1. Cyanate ion, OCN" 

The 2165 cm"' (4.62 jjm) peak, commonly referred to the 
"XCN" band, was first detected toward the massive protostar 
W33A (Soifer et al. 1979). Its presence was the first observa- 
tional indication that complex chemistry (triggered by UV/X- 
ray photons or ion/electron bombardments) could be occurring 
in the interstellar ice mantles. To date this feature has already 
been observed in more than 30 deeply embedded mostly high- 
mass young stellar objects (Tegler et al. 1993; Demyk et al. 
1998; Pendleton et al. 1999; Gibb et al. 2000; Keane et al. 2001; 
van Broekhuizen et al. 2005; Pontoppidan et al. 2003; Whittet 
et al. 2001 ; Gibb et al. 2004) and several galactic center sources 
(Chiar et al. 2002; Spoon et al. 2003). 

The "XCN" band has been studied extensively in the labo- 
ratory, being easily produced by proton irradiation (e.g. Moore 
et al. 1983; Hudson et al. 2001), vacuum ultraviolet pho- 
tolysis (e.g. van Broekhuizen et al. 2005; Lacy et al. 1984; 
Grim & Greenberb 1987; Demyk et al. 1998), or thermal an- 
nealing of interstellar ices analogues (Raunier et al. 2003; 
van Broekhuizen et al. 2004). Along the years, several can- 
didates have been proposed to justify the XCN band, includ- 
ing some isonitriles (X-NsC) and cyanates (X-0-C=N) (see 
Pendleton et al. 1999 for review). However, laboratory ex- 
periments involving isotopic substitution proved unequivocally 
that the responsible for the 2165 cm"' feature is the cyanate 
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Fig. 5. Comparison between IR spectra of interstellar and lab- 
oratory ices. The highest curve is the infrared spectra of proto- 
stellar source W33A obtained by Infrared Space Observatory 
(ISO). Lower traces indicate lab spectra of H20:NH3:CO ices 
after processing by: (a) UV photons (Hudson & Moore 2000); 
(b) 0.8 MeV protons (Hudson & Moore 2000) and (c) 46 MeV 
Ni ions (this work). 

ion, OCN^ (Schutte & Greenberg 1997; Bernstein et al. 2000; 
Novozamsky et al. 2001; Palumbo et al. 2000). 

The present work is the first report on the OCN" formation 
by heavy, highly charged and energetic ions (46 MeV ^**Ni'^^), 
as an attempt to simulate the chemistry and the physical chem- 
istry induced by cosmic rays inside dense regions of interstellar 
medium like dense molecular clouds or protoplanetary disks. 

Following Hudson et al. (2001) a possible route for the pro- 
duction of OCN , from the radiolysis of interstellar ice ana- 
logues, occurs via the reaction of NH or NH2 radical with CO, 
leading to the formation of isocyanic acid molecule HNCO: 

NH3 (NH^ + e") or (NH2 + H) (8) 

NH3 + NHJ ^ NH| + NH2 (9) 

NH2 -H CO ^ H -H HNCO or NH + CO HNCO (10) 

The reaction between isocyanic acid and ammonia (via 
acid-base reaction) produces ammonium, NH^, and OCN : 

HNCO + NH3 ^ NH+ + OCN (11) 

In the infrared spectra of the irradiated ices by heavy ions, 
beyond the OCN" band, the HNCO and NH| species are also 
observed at frequencies around 2265 and 1490 cm"\ respec- 
tively. 

A comparison between the infrared spectra of interstel- 
lar ices around the embedded protostar W33A with irradi- 
ated ammonia-water-CO mixtures from laboratory is given in 
Fig. |5] The highest curve presents the IR spectrum recorded 
by Infrared Space Observatory (ISO) toward the protostel- 
lar object W33A. The data have been extracted from ISO 
database]^. Curve (a) shows the IR spectrum of H20:NH3:CO 

^ http://iso.esac.esa.int/ida 



(1:0.2:0.2) ice at 15 K processed 10 minutes of UV pho- 
tons from hydrogen-discharge lamp (Hudson & Moore 2000). 
Curve (b) presents the IR spectrum of H20:NH3:CO (1:0.2:0.2) 
ice at 15 K irradiated by 0.8 MeV protons at a dose of 19 
eV molec ' (Hudson & Moore 2000). The lowest curve (c) 
shows the infrared spectrum of H20:NH3:CO (1:0.6:0.4) ice 
at 13 K observed after the bombardment by 46 MeV Ni ions at 
fluences around 2x10'^ ions cm"^ (this work). The radiolysis 
of H20:NH3:CO performed by the present study, using heavy 
and energetic cosmic rays, reproduces very well the OCN" 
(2165 cm"') and CO (2139 cm"') infrared features observed 
in the W33A spectrum. Other spectral features like the broad 
IR peaks observed at 1650 and 1450 cm"' are also similar to 
the astronomical source. 

4.3.2. Ammonium, NH+ 

The broad V4 vibration mode of NHJ around 1490 cm"' is also 
observed after the heavy ion bombardment on both ammonia- 
containing ices (see Figs |2^-b). This line also has been ob- 
served from similar irradiation experiments involving photol- 
ysis and radiolysis of ammonia-containing ices (e.g. Demyk et 
al. 1998; Hudson et al. 2001; Hudson & Moore 2000). As dis- 
cussed by Demyk et al. (1998), the band strength of the NHJ 
V4 line is very uncertain (due to the broad profile in the as- 
trophysical ice analogs). This fact, combined with a possible 
contamination from other infrared lines, makes this particular 
column density determination very imprecise. 

The NHJ formation directly from the processing of ammo- 
nia is given by Eqs.[8]and|9] The reaction pathways due to the 
processing of H20:NH3:CO ice are given by Eqs. [TOlandfTTI 
As discussed by Moore et al. (2007), the presence of water in 
interstellar ices also contributes to convert NH3 into NH^ : 

H2O H2O+ H- e" (12) 

H2O+ + H2O ^ OH -H H3O+ (13) 

H2O+ H- NH3 ^ OH -H NH| (14) 

H3O+ + NH3 ^ H2O + NH^ (15) 

Besides the observation of the OCN" ion, the occurrence 
of ammonium at several protostellar disks was also observed 
(Gibb et al 2000; Boogert et al. 2004; Gibb et al. 2004; Demyk 
et al. 1998 and references therein) and possibilly for several so- 
lar system ices (Moore et al. 2007). Following Muiioz Caro & 
Shutte (2003), NH^ is one of the species involved in the forma- 
tion of highly complex molecules like hexamethylenetetramine 
HMT, (CH2)6N4, which was observed in experiments involv- 
ing UV photoprocessing of interstellar ammonia-containing ice 
analogs. 

4.3.3. Other species 

The infrared spectra of H20:NH3:CO ice (1:0.6:0.4) from 2400 
to 1200 cm"' acquired during the slow sample heating, from 
13 K up to room temperature, are shown in Fig. |6^. For each 
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spectrum the sample temperature is indicated. The spectra have 
an offset of 0.02 for better visualization. The CO line at 2139 
cm"' decreases drastically in a temperature between 13 and 52 
K. At 150 K the N2O peak (~ 2234 cm ') increases sHghtly. 
The unidentified feature around 1450 cm ' is not seen any- 
more above this temperature. At 200 K all the water, ammonia 
and CO2 ices have been sublimated and their peaks have dis- 
appeared. Although very attenuated, the OCN" peak (~2169 
cm"') is still present at 200 K. At this temperatures and higher, 
we observe the appearance of new unidentified line around 
2147 cm"' (between the OCN" and CO peaks). This peak is 
tentatively assigned to aliphatic isocyanide R-N=C molecules 
(2150 cm"') as observed by Imanaka et al. (2004) in UV pho- 
tolysis of N2:CH4 ices. 

A comparison between the irradiated ice at 1 3 K and the or- 
ganic residue at 300 K is shown at Fig.|6j3. Tentative band attri- 
bution is indicated. At room temperature, a large structure (pos- 
sibly two overlapping peaks) is observed around 2210 cm"'. 
Since these features are close to C-N stretch vibration band, 
we suggest that they may be due to organic residues contain- 
ing CN bearing molecules. Vertical dashed lines indicate the 
possible frequency of some vibration modes of the zwitterionic 
glycine (NH+CH2COO"). 

The observed infrared bands produced by the radiolysis of 
the H20:NH3:CO ice (1:0.6:0.4) at 13 K by heavy, charged 
and energetic ions in the 2400 to 1200 cm"' range are listed 
in TablelD The lines observed after the sample heating to room 
temperature are also shown. The assignments are proposed by 
comparing the spectra with those of similar experiments by 
photolysis or radiolysis of interstellar ice analogs (Hudson et 
al. 2001; Jamieson et al. 2005; Hudson & Moore 2000; van 
Broekhuizen et al. 2005; Munoz Caro and Shutte 2003; Holtom 
et al. 2005; Demyk et al. 1998; Imanaka et al. 2004). The 
low frequency lines are tentatively assigned to N vO, molecules 
(Jamieson et al. 2005). 

At room temperature, the peaks around 1593, 1506, 1415, 
1338 cm"' may possibly be assigned to different vibration 
modes of zwitterionic glycine (NH:|CH2COO"). This amino 
acid was observed among the residues coming from 5 keV elec- 
tron bombardment of CH3NH2:C02 (30:1) ice (Holtom et al. 
2005) at 10 K. Zwiterionic glycine has been also observed af- 
ter iiTadiation of NH3:CH3COOD (1:1) ice by low energy elec- 
trons (Lafosse et al. 2006) at 25 K. The authors have observed 
zwitterionic glycine even at 25 K, without any subsequent ther- 
mal activation. 

The presence of amino acid precursors has been observed 
among the organic residues due to photolysis and radiolysis 
of several ammonia-containing interstellar ice analogs. The in- 
vestigation employed vacuum ultraviolet photons from syn- 
chrotron light sources (Chen et al. 2008; Nuevo et al. 2007; 
Nuevo et al. 2006), hydrogen UV-lamp (Elsila et al. 2007; 
Munoz-Caro and Schuttle 2003; Bernstein et al. 2002; Munoz 
caro et al. 2002) and recently, by ion bombardment (Takano et 
al. 2007; Kobayashi et al. 2008). 

The IR peak centered around 1370 cm"' observed in 
both 13 K and 300 K spectra was tentatively assigned 
to hexamethylenetetramine - HMT, (CH2)6N4, since a sim- 
ilar feature was also observed in the UV photolysis of 
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Fig. 6. (a) Infrared spectra of H20:NH3:CO ice (1:0.6:0.4) from 
2400 to 1200 cm"' during heating up to room temperature. The 
sample temperature for each spectrum is given. Each spectrum 
has an offset of 0.02 for better visualization, (b) Comparison 
between the irradiated ice at 13 K (top spectrum) and the 300 
K residue (bottom spectrum). Vertical dashed lines indicate the 
frequencies of some vibration modes of zwitterionic glycine 
(NH+CH2COO ). 

H20:NH3:CH30H:CO:C02 ice by Munoz Caro and Shutte 
(2003). The observation of NH^ in the same organic residue 
reinforces the evidence for the presence of HMT, since the am- 
monium is an essential component for its formation. 

The presence of these complex molecules in the organic 
residue from the bombardment of interstellar ice analogs by en- 
ergetic and heavy ions has to be confirmed. However, the result 
of this experiment suggests that even deeply inside molecular 
cores or other denser regions like protoplanetary disks, the in- 
terstellar grains are being highly transformed. The produced or- 
ganic molecules, trapped into and onto dust grains, meteoroids 
and comets, could be delivered into the planets/moons possi- 
bly allowing pre-biotic chemistry in such environments where 
water is also found in liquid state. 

4.4. Astrophysical implication 

Inside dense regions of interstellar medium, where direct stel- 
lar UV is extremely attenuated. X-rays (despite some attenua- 
tion degree) and cosmic rays are the main drivers for the gas 



Pilling et al.: Radiolysis of ammonia-containing ices by heavy and energetic cosmic rays 



9 



Typical dense cloud grain 



Hea^y cosmic ray 



Confined sputtered 
molecules 



Ion track 



Sputtered 
molecules 



I 




I Time scale \ 1 e-^' 

10 s < I e-^^ 
I (cadiA) ^ 



Ion track 



\.im 



100 nm 
lOnm 



Re-adsorbed 
molecules ' 



Water-rich ice mantle 
Silicate/Carbon grain 




Sputtered 
light species 



Sputtered t h,o 
I molecular "2° ' f oh- * co* f 

I clusters --^f . _ f • \ ^ | 



Sputtered 
species 



Fig. 7. Schematic view showing the interaction between heavy ion cosmic ray with a typical interstellar grain inside dense 
clouds. The ion track along the coagulate sub-micron size grains, the grain mantles, the processed and the sputtered molecules 
are indicated. Figure insets were adapted from Andrade et al. (2008) and indicate the physical-chemical changes on the grain 
mantle due to the impact with heavy ion. See details in text. 



phase and grain surface chemistry. For example, the penetra- 
tion depth of 46 MeV Ni ions and equivelocity protons inside 
water ice is about 1 and 3 orders of magnitude higher than for 
500 and 50 eV photons (considering the photoabsorption cross 
section from Chan et al. (1993) and McLaren et al. (1987), re- 
spectively). In these environments the dust grains can reach a 
size of the order of a micron, as a coagulation of sub-micron 
grains with ice-rich mantles of tens of nanometers (Mathis 
1990, Mathis, Rumpl, Nordsieck 1977; Li & Greenberg 1997; 
Chokshi, Tielens & Hollenbach 1993). However, the presence 
of larger grains (radii within one to tens of microns) has also 
been suggested (Witt, Smith & Dwek 2001 ; Taylor Baggaley & 
Steel 1996;Griinetal. 1994). 

FigureQillustrates, as a schematic view, the interaction be- 
tween heavy ion cosmic rays with a typical interstellar grain 
inside dense clouds. The ion track (along the coagulated sub- 
micron size grains), the grain mantles, the processed material 
and the sputtered molecules are indicated. The inset of the fig- 
ure was adapted from Andrade et al. (2008) and shows the dif- 
ferents stages (times scales) of physico-chemical changes on 
the grain mantel caused by the heavy ion impact. At the per- 
tinent heavy ion velocities, the projectile traverses atomic dis- 
tances of the grain mantle in 10"'^ s. The energetic cosmic ray 
analog ionizes and excites molecules along its trajectory inside 
the solid (inset fig.|7^). The secondary electrons (6) move away 
from the projectile trajectory, generating a positive infratrack in 
the center and a negative ultratrack around. Afterwards, molec- 
ular dissociations and chemical reactions occur, causing sput- 
tering of light ions such as (inset fig.|7];). Following Iza et 
al. (2006) the sputtering occurs mainly due to long lived re- 
pulsive electronic excitation which lead to atomic or molecular 
motion. 



After the Coulomb repulsion of ionizing species, the next 
stage is the track relaxation, when preformed or newly formed 
larger chemical species (molecular clusters) are desorbed (in- 
set fig.|7}l) at timescales around ~ 10"'° s. Shock waves trans- 
fer momentum to material at the ice surface which desorb as 
molecules from a site located as far as 3-5 nm of the ion track. 

Due to the high porosity of grains inside dense clouds 
(coagulated from sub-micron heterogenic grains), some sput- 
tered species inside the grain pores become confined, being 
re-adsorbed on the grains at a later time. The re-adsorption of 
sputtered species (processed or not) my also occur at the outer 
layer of the micron size grain mantle. The projectile energy de- 
posited along the interfaces between the water-rich mantle and 
the silicate/carbon grain may promote chemical mixing within 
these two different chemical regions. 

Inside dense interstellar regions, the 1 MeV proton flux is 
around 1 proton cm"^ s ' (Mennella et al. 2003; Morfill et al. 
1976). Inside these regions the iron to proton ratio is about 
1.6 xlO and the iron fraction of the canonical cosmic ray 
flux is estimated to be 2 x 10 cm"^ s ' (Roberts et al. 2007; 
Leger et al. 1985). Following Drury et al. (1999), heavy ions 
with 12 < Z < 29 (e.g. from Mg to Cu) present roughly the 
same abundance, while those with Z > 30 are more than 3 or- 
ders of magnitude less abundant. The average cosmic ray flux 
composed by heavy nuclei ions (Z > 12) is estimated to be 
<pHCR ~ 5 X 10~^ cm"^ s"'. Despite this small value, the sput- 
tering promoted by heavy and energetic ions can be more than 
4 orders of magnitude higher than ones promoted by the abun- 
dant low energy protons (Leger et al. (1985), Seperuelo Duarte 
et al. (2009a,b)). Moreover, in regions where the ionizing stel- 
lar photons are very attenuated, the sputtering promoted by cos- 
mic rays becomes a competitive desorption process. This issue 
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Table 4. Assignment of infrared absorption features produced 
by the radiolysis of the H20:NH3:CO ice (1:0.6:0.4) by 46 
MeV 58Nii3+ iojis at 13 K and after warming to 300 K. 



Frequency 


Wavelength 


Temp. 


Molecule 


Notes 


(cm-') 


(fim) 


(K) 






2233 


4.48 


13 


NjO 


[1,2] 


2218-2200 


4.51-4.54 


300 


nitrilest 


[8] 


2168 


4.61 


13, 300 


OCN- 


[1,3,4,7] 


2147 


4.66 


300 


aliph. isocyanidet 


[8] 


~ 2112 


4.73 


300 


NCOj 


[2] 


1725 


5.80 


300 


ester"'" 


[5] 


1683 


5.94 


300 


amides^ 


[5] 


1652 


6.05 


300 


asym-NiO^ 


[2] 


1637 


6.11 


13 


7 




1593 


6.28 


300 


NH+CHjCOO-t 


[6] 


1558 


6.42 


300 


? 




1533 


6.52 


300 


? 




1506 


6.64 


300 


NH+CH2COO-' 


[6] 


-1490 


6.71 


13 


NH+ 


[1,3,7] 


1474 


6.78 


13 


NO^ 


[2] 


1440 


6.94 


13 


NH+CHjCOO-t 


[6] 


1415 


7.07 


300 


NH+CH2COO-' 


[6] 


~ 1370 


7.30 


13, 300 


HMT+ 


[5] 








HCOO- 


[5,9] 


-1338 


7.47 


13, 300 


NH+CH.COO-t 


[6] 








NHjCHzCOQ-^ 


[6] 








HCOO 


[9] 


1305 


7.66 


13 


N205;N20i 


[2] 


1283 


7.80 


300 


N2OI" 


[2] 



'Tentative assignment. 

[1] Proton bombardment of several ices (Hudson et al. 2001); [2] 
Electron bombardment of N2:C02 (Jamieson et al. 2005); [3] Hudson 
& Moore 2000; [4] van Broekhuizen et al. 2005; [5] UV photolysis 
of H20:NH3:CH30H:CO:C02 ice (Munoz Caro and Shutte 2003); 
[6] Electron bombardment of CH3NH2:C02 ice (Holtom et al. 2005); 
[7] UV photolysis of ammonia-containing ices (Demyk et al. 1998); 
[8] UV photolysis of N2:CH4 ices at various pressures (Imanaka et 
al. 2004); [9] proton bombardment of H20:NH3:CO ice (Hudson & 
Moore 2001). 

Table 5. Average dissociation rate (R) and the corresponding 
half life (fi/2) of water, ammonia and CO molecules under 
bombardment by 46 MeV Ni ions; similar results are expected 
for heavy cosmic rays interaction with a interstellar ammonia- 
containing (H20:NH3:CO (1:0.6:0.4)) grain. 



Species 


Relative initial 


R 


tl/2 




abundance 


(10"'^ s-') 


{10^ years) 


H2O 


50% 


8 


2 


NH3 


30% 


7 


3 


CO 


20% 


9 


2 



should be taken in consideration in future chemical models of 
inner regions of molecular clouds where a considerable fraction 
of excited molecules are released into the gas phase by heavy 
ion impact. 



A rough estimation of half lives of molecular species in 
space can be carried out by considering only molecular disso- 
ciation effects of cosmic rays. The average dissociation rate, R, 
of a given species /, induced by the interstellar heavy cosmic 
ray field is approximative ly given by: 

Ri ~ <pHCR X o-d.i [s^'l (16) 

where cr,/ ,■ is the average dissociation cross section of a frozen 
species / due to heavy cosmic rays and (pncR is the average 
heavy ion cosmic ray flux (ions cm^^i"'). The half-life, fi/2, 
for a given species / may be obtained from Eq.[T6]by writing: 

In 2 

fi/2,, = — [s] (17) 

Ki 

which does not depend on the molecular number density in the 
ice. The dissociation cross section depends on the projectile 
energy: for Ni ions, crj j is typically 2x10"'-' cm^ at 46 MeV, 
decreasing to 0.3 xlO"'"' cm^ at 537 MeV (Seperuelo Duarte 
et al. 2009b). Nevertheless, assuming the values of o-^ ,- for Ni 
ions at 46 MeV as the average value and <pHCR ~ 5 x 10"^ ions 
cm"^ s"' as the heavy ion cosmic ray flux, R can be estimated. 

Table |5]presents the average dissociation rate and half life 
of water, ammonia and CO molecules bombarded by heavy 
ions The estimated half life for the studied species over a typ- 
ical heavy cosmic ray flux is about 2-3 xlO^ years. This value 
is in a good agreement with the half lives of typical molecu- 
lar clouds and protostellar clouds (e.g. Millar & Nejad 1985; 
Takeuchi et al. 2005). Better estimation for molecular dissocia- 
tion rates and half lives promoted by the direct impact of heavy 
cosmic rays on interstellar ices depends on more accurate mea- 
surements/determination of the averaged heavy ion cosmic ray 
flux inside these dense environments. 

5. Summary and conclusions 

Experimental study on the interaction of heavy, highly charged 
and energetic ions (46 MeV ^^Ni'-'^) with ammonia-containing 
ices H20:NH3 (1:0.5) and H20:NH3:CO (1:0.6:0.4) is per- 
formed to simulate the physical chemistry induced by cosmic 
rays inside dense regions of interstellar medium like dense 
molecular clouds or protoplanetary disks. Our conclusions are: 

1. At the beginning of the iiTadiation, the ices become com- 
pact as revealed by the particular dependence of the column 
density on the fluence. Measuring the degree of compaction 
through the integrated area of the OH dangling bond feature 
(~ 3650 cm"') from of water molecules trapped into ice mi- 
cropores, the eff'ect promoted by the impact of heavy ions 
seems to be at least 3 orders of magnitude higher than the 
one promoted by (0.8 MeV) protons. 

2. The infrared spectra of the irradiated samples present the 
production of several new species, including OCN , HNCO 
and NH+. The OCN" band (2169 cm"') is observed even 
at room temperature. During the ice heating a broad fea- 
ture was observed around 2218-2200 cm"' and was tenta- 
tively assigned to non-volatile nitriles. At room tempera- 
ture an intense and sharp feature was observed at (~ 2150 
cm"') and was tentatively assigned to aliphatic isocyanide 
(-2150 cm"'). 
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3. The spectra of the irradiated HjOiNHjiCO ice (1:0.6:0.4) 
at room temperature reveal several bands which are tenta- 
tively assigned to zwitterionic glycine and HMT. Despite 
the presence of these complex molecules in the organic 
residue from the bombardments of interstellar ice analogs 
by energetic and heavy ions has to be confirmed. 

4. The obtained value for the dissociation cross section of wa- 
ter, ammonia and carbon monoxide (in the tertiary mix- 
ture ice) due to heavy cosmic ray analogs are ~ 2 x 10 
1.4x10"^^ and 1.9x10"'^ cm^, respectively. These values 
seem not to be affected by small (<30%) changes in the 
initial relative molecular abundances. 

5. In the presence of a typical heavy cosmic ray field the es- 
timated half lives for the studied species is about 2-3 xlO'' 
years. This value is in a good agreement with the half fives 
of typical molecular clouds and protosteUar clouds. 

Although they represent only a small fraction (~ 1%) of 
the cosmic rays flux, some effects (e.g. molecular sputtering 
and ice compaction) promoted by heavy ions on interstellar ice 
grains are more intense than those promoted by protons. This 
should be taken under consideration in future chemical models 
of inner regions of molecular clouds since heavy ions wiU con- 
tribute molecules to the gas phase and trigger chemical surface 
and bulk reactions. 
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